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Abstract 

Unlike convex case, a local equilibrium point of a nonconvex Nash-Cournot oligopolistic 
, equilibrium problem may not be a global one. Finding such a local equilibrium pomt or even 

a stationary point of this problem is not an easy task. This paper deals with a numerical 
method for Nash-Cournot equilibrium models involving nonconvex cost functions. We develop 
■ a local method to compute a stationary point of this class of problems. The convergence of 

the algorithm is proved and its complexity is estimated under certain assumptions. Numerical 
examples are implemented to illustrate the convergence behavior of the proposed algorithm. 
Keywords. Nonconvex Cournot-Nash models, splitting proximal point method, local equilib- 
ria, gradient mapping. 
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1 Introduction 



Nash-Cournot oligopolistic equilibrium models have been widely applied in economics, electricity 
markets, transportation, networks as well as in environments. Such a model can be formulated as 
a game model, where each player has a profit function which can be expressed as the difference of 
^ I its price and cost functions. In classical models, the price function is affine while the cost function 

' is assumed to be convex. In this case, a local equilibrium point is also a global one. Mathematical 

programming and variational inequality approaches can be used to treat this problem (see, e.g., 

[anisiiniiii]). 

In practical models, the cost per a unit usually decreases as the production level increases. This 
1^ ' situation requires a nonconvex function to represent the production cost of the model. In [TS], a 

global optimization algorithm has been developed to find a global equilibrium point to the Nash- 
Cournot oligopolistic equilibrium market model involving piecewise concave cost functions. However, 
global algorithms only work well for the problems of moderate size, while it becomes intractable when 
the size of the problem increases, except for special structures are exploited. 
I In this paper, we continue the work in [15] by proposing a local solution method for finding 

H ' a stationary point of Nash-Cournot equilibrium models with nonconvex (not necessarily concave) 

. . .' cost functions. We consider a Nash-Cournot model involving an affine price function and nonconvex 

smooth production cost functions. With this structure, the cost function of the model can be 
decomposed as the sum of a convex quadratic function and a nonconvex smooth function. Then, 
we develop a local method for finding a stationary point of such a model. The method is called a 
splitting proximal point algorithm. The main idea of this algorithm is to preserve the convexity of 
the problem while convexifies the nonconvex part by linearizing it around each iteration point. 

Proximal point methods have been well developed in optimization as well as in nonlinear analysis. 
Myriad of research papers concerned to these methods were published (see, e.g., [TTJ [T^l [221 125] 
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and the references quoted therein). However, all these papers only deal with a class of monotone 
problems. Recently, Pennanen in [21J extended the proximal point method to nonmonotone cases 
for solving variational inequality problems. Lewis further generalized this algorithm in a unified 
framework using the prox- regularity concept f23^ • A main point of the proximal point methods is 
to choose the proximal parameter sequence. This affects to the performance of the algorithm as 
well as its global convergence behavior. Giiler in [5] investigated the rate of global convergence of 
the classical proximal point methods for convex programs. The worst case complexity bound of this 
method is 0(l/fc), where k the number of iterations. The author further accelerated the classical 
proximal point method to get a better complexity bound for the convex programming problems, 
precisely, 0{l/k^) [7] by using the idea of Nesterov for gradient methods [18]. A splitting proximal 
point are also developed in many research papers. Such a method was applied to optimization 
problems by Mine and Fukushima in |13| and, recently, by Nesterov in [17] . 

This paper contributes a new local method for finding a stationary point of a Nash-Cournot 
equilibrium models involving the nonconvex cost functions. The algorithm called splitting proximal 
point method is presented and its convergence is investigated. The global worst case complexity 
bound is also provided, which is 0{1/Vk), where k is the iteration counter. To our knowledge, this 
is the first estimation proposed to Nash-Cournot equilibrium model. 

The rest of the paper is organized as follows. Section [2] presents a formulation of a Nash-Cournot 
oligopolistic equilibrium model involving nonconvex cost function. The problem is reformulated as 
a mixed variational inequality. In Section [3j we define three concepts including local and global 
equilibria, and stationary points of the Nash-Cournot equilibrium model. Section |4] deals with a 
gradient mapping and its properties. The splitting algorithm is described in Section [SJ where its 
convergence is proved and the worst-case complexity is estimated. Two numerical examples are 
implemented in the last section. 

2 Mixed variational inequality formulation 

We consider a Nash-Cournot oligopolistic equilibrium market models with n-firms producing a com- 
mon homogeneous commodity in a non-cooperative fashion. The price p of production depends on 
the total quantity a :— X]"=i of the commodity. Let hi(xi) denote the cost of the firm i when its 
production level is Xi. Suppose that the profit of firm i is given as 

n 

fi{xi,- ■ ■ ,Xn) ■■= Xip(^Xi) - hi{xi), i = l,...,n, (2.1) 

1=1 

where hi is the cost function of the firm i which is assumed to only depend on its production level. 

Let Ci C M (i = 1, . . . , n) denote the strategy set of the firm i, which is assumed to be closed and 
convex. Each firm seeks on its strategy set to maximize the profit by choosing the corresponding 
production level under the presumption that the production of the other firms are parametric inputs. 
In this context, a Nash equilibrium is a production pattern in which no firm can increase its profit by 
changing its controlled variables. Thus under this equilibrium concept, each firm determines its best 
response given other firms' actions. Mathematically, a point x* = {xl, . . . , x* )"^ G C := Ci x • • • x C„ 
is said to be a Nash equilibrium point if 

Mx*i,-- ■ ,a;*_i,2/.,,a;*+i, • • • ,x*J < f,{xl, ■ ■ ■ ,x*J, \/yi e Ci (i = 1, . . . ,n). (2.2) 

When hi is affine, this market problem can be formulated as a special Nash equilibrium problem 
in the n-person non-cooperative game model, which in turns is a strongly monotone variational 
inequality (see, e.g., [9J). 
Let us define 

n 

*(a;,2/) := - ^Mxi,- ■ ■ ,Xi-i,y^,Xi+i, ■ ■ ■ ,x„), (2.3) 

i=l 

and 

0(x,?/) :=*(x,2/)-*(a;,a;). (2.4) 
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Then, as proven in |5] , the problem of finding an equilibrium point of this model can be reformulated 
as the following equilibrium problem: 

Find X* &C such that: (j){x\y) > for all y £ C. (EP) 

This generalized setting was proposed by Blum and Oettli in [5] (see also in [TFl). 

In classical Cournot models [U [SJ , the price and the cost functions for each firm are assumed to 
be affine and given as follows: 

n 

p{a) — ao — (3a, ao > 0, /? > 0, with a — "^^Xi, (2-5) 

1=1 

h,{xi) ^ fi^Xi + Ml > 0, 6>0 = (2.6) 
In this case, using (pH . ([^ . (^3)) and ((^ . it is easy to check that 

Hx, v) = {Bx + n-af{y-x) + ^y'^ By - ^x'^Bx, 

where 

/3 P ■■■ P' 
P p ■■■ p 

1 

p p p ■■■ 

a = (ao, . . . ,ao)^, and = (^i, . . . , ^„)^. 

Then the problem of finding a Nash equilibrium point can be formulated as a mixed variational 
inequality of the form: 

Find X* £ C such that: 
{Bx*+^i-af{y-x*) + ^y'^By^^{x*fBx* > 0, Vy G C. (2.7) 

Let Q := B + B . Since /? > and matrices i? and B are symmetric, it is clear that Q is symmetric 
and positive definite. This mixed variational inequality can be reformulated equivalently to the 
following strongly convex quadratic programming problem: 

min {ix^Qa; + (m - a)'^x}. (QP) 

x£U 2 

Hence, problem ( |QP[ ) has a unique optimal solution, which is also the unique equilibrium point of 
the classical oligopolistic equilibrium market model. 

The oligopolistic market equilibrium models, where the profit functions (« = 1, . . . , n) of each 
firm is assumed to be differentiable and convex with respect to its production level Xi while the 
other production levels are fixed, are studied in |3| (see also [5]). This convex model is reformulated 
equivalently to a monotone variational inequality. 

In practical models, the production cost function hi assumed to be affine is no longer satisfied. 
Since the cost per a unit of the action does decrease when the quantity of the commodity exceeds 
a certain amount. Taking into account this fact, in the sequel, we consider market equilibrium 
models where the cost function h may not be convex, whereas the price function is affine as in (|2.5p . 
Typically, the cost function h is given as: 

n 

h{x) ■.= Y,'^'^i^^)> (2-8) 

1=1 

where hi {i = 1, ... ,n) is differentiable and nonconvex. 
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Let us denote by 



F(x) := Bx — a 
(p(x) := g{x) - h{x), 

where a := a — /i, g{x) :— ^x'^ Bx and h[x) defined as (j2.8p . 

Obviously, matrix B is symmetric, using the same notation i?, B and a as in (j2.7p . we can 
formulate the nonconvex Nash-Cournot equilibrium model as a mixed variational inequality: 

Find X* £ C such that: 

F{x*f{y - X*) + ip{y) - ip{x*) > for all y e C. (ncMVIP) 

Note that mixed variational inequality problems of the form (jncMVIPp , where (f is convex, i.e. h is 
concave, were extensively studied in the literature (see, e.g., pj |4l 13 [9l [TOl [IH [191 124]). 

Remark 2.1. If the function (p is convex and differentiable then problem (|ncMVIP[) can be reformu- 
lated equivalently to a classical variational inequality problem. More generally, it can be converted 
to a generalized variational inequality problem when ip is convex and subdifferentiable (see, e.g. 
[9]). However, the mixed variational inequality (jncMVIPl) can not equivalently transform into a 
variational inequality if ip is nonconvex. 

Note, if we define (j){x, y) := F{x)^{y — x) + p{y) — p{x) then problem (jncMVIPp coincides with 
a nonconvex equilibrium problem of the form (jEPp . 

3 Local equilibria and stationary points 

Unlike the convex case, if the cost function ip of the problem (jncMVIPp is nonconvex, it may not 
have a global equilibria even if C is compact, and F and p are continuous. Indeed, let us consider 
C :— [—1, 1] C K, F{x) :— x and ip{x) = —\x^, which is concave, then F{x)'^{y — x) + (p{y) — (p{x) — 
~^{y — x)'^. Therefore, problem (jncMVIPP corresponding to this function has no solution. 

For a given x G C, let M{x, r) be an open ball centered at x of radius r > in M". Borrowing 
the concepts from classical optimization, we firstly propose a local equilibria and critical points (or 
stationary points) of the mixed variational inequality (IncMVIPp . 

Definition 3.1. A point x* G C is called a local solution (or local equilibria) to (jncMVIPl) if there 
exists a ball M(x*,r) such that 

F{x*f{y~x*) + p{y) - ip{x*) > for all y € CnB(a;*,r). (3.1) 

If C C M{x*,r) then x* is called a global solution (or global equilibria) to (jncMVIPp . 

Let 

^{x; y) := {Bx ~ af{y - x) + ^y'^By - h{y). (3.2) 

We consider a function m : C x R++ M and a mapping S : C x R++ =^ 2'-' defined as follows: 

m(x; r) := min {^{x;y) \ y e C nM{x,r)} , (3.3) 
S{x; r) := argmin {?A(a:; y) | y G C n B(x, r)} , (3.4) 

where B(a;, r) stands for the closure of the open ball M{x, r). As usual, we refer to m as a local gap 
function for problem (jncMVIPp . Obviously, if h is continuous and C is compact then the function 
m as well as the mapping S are well-defined. If h is concave then S is reduced to a single valued 
mapping due to the symmetric and positive definiteness of B. 

The following proposition gives a necessary and sufficient conditions for a point to be a local or 
global solution to (|ncMVIPp . 
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Proposition 3.1. The function m defined by (j3.3l) satisfies m{x) < for all x £ C . Moreover, the 
following statements are equivalent: 

a) X* is a local solution to (|ncMVIP[) .' 

h) There exists f > and x* G C such that m(x*; f) = 0; 

c) There exists r > and x* G C such that x* G S{x*; f). 

Proof. Note that ip{x,x) = for all a; G C, and for any x G C and f > 0, x G CD M{x,f). 
Therefore, from the definitfon (j3.3l) of to, it is clear that, with a given r > 0, TO(a;; f) = 
min {■;/'( x, y) | y G C nB(a;, f)} < ip{x,x) = for every cc G C. The equivalence bewteen b) and 
c) is trivial. We only prove that a) is equivalent to b). 

Suppose that there exists f > and x* G C such that m{x* , r) ~ 0. It follows from the definition 
of TO that Tp{x*,y) > 'ip{x*,x*) = for all x G B(x*,f) n C. In particular, ip{x*,y) > for all 
y E C O M{x* , f). Thus x* is a local equilibria of (jncMVIPl) . Conversely, if x* is a local equilibria of 
(jncMVIPp then there exists a neighbourhood 8(2:*, r) such that r > and 

F{x*f{y-x*) + ip{y) - ip{x*) > F{x*f{x* - x*) + ip{x*) - ip{x*) = 0, Vy G Cn'B{x*,r). 

Since r > and M{x*,r) C M", there exists < f such that l(x*,f) C B{x*,r). This relation 
shows that the last inequality holds for all y G C n R{x* ,f). Therefore, m{x* , f) — 0. □ 

Clearly, if the conclusions of Proposition 13. II hold for a given r > and C C M{x* ,f) then x* is 
a global solution to (jncMVIPp . 
Next, let us define 

Tc{x) {d ■.= t{y-x) \y€C, t>0}, (3.5) 

a cone of all feasible directions of C starting from x G C. The dual cone of Tc{x) is the normal 
cone of C at x which is defined as 

iV<.(a;) •= /{^^^" I (36) 
1 0, otherwise. 

By Proposition l3.11 a point x G C is a local solution to (jncMVIPp if and only if it solves the following 
optimization problem: 

X G argmin|(i?x - a)'^{y - x) + ^y'^ By - h{y) \ y G Cnl(x,f)| , (3.7) 

for some f > 0. Since h is not necessarily concave, finding such a point x satisfying (13. 7|) . in general, 
is not an easy task. In this paper, we concentrate in finding a stationary point rather than local 
equilibrium. We develop a method to find such a point for (|ncMVIP|) . Borrowing the concept of 
stationary points in optimization, we define a stationary point (or a critical point) for the mixed 
variational inequality (jncMVIPp as follows. 

Definition 3.2. A point x G C is called a stationary point (or critical point) to the problem 
(|ncMVIP|) if 

e Qx~ a-Vhix) + Ncix), (3.8) 
where Nc is defined by (13.61) and Q -.^ B + B. 

Since Nc is a cone, for any c > 0, the inclusion p.Sp is equivalent to 

Oec[{B + B)x-a-Vh{x)]+Nc{x). (3.9) 

Let 

D(j){x; d) := [{B + B)x - & - Vh{x)fd, (3.10) 
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for any x € C and d G J-c{x). Then the condition (I3.8P is equivalent to 

D4>ix*;d) >0, ydeTcix*). (3.11) 

Let us denote by S* the set of stationary points of (jncMVIPp . The fohowing lemma shows that 
every local equilibria of problem (jncMVIPl) is its stationary point. The proof is simple and short, 
we present here for reading convenience. 

Lemma 3.1. Suppose that h is continuously dijjerentiable on its domain. Then, every local equilibria 
is a stationary point of problem ()ncMVIPp . 

Proof. Suppose that x* is a local equilibria of (|ncMVIP|) . Then there exists a neighborhood M{x* ,r) 
of X* such that F{x*)'^{y — x*) + ip{y) — 'f{x*) > for all y € M{x* ,r) n C. According to Proposition 
13.11 this requirement is equivalent to x* is a local solution of 

minl^iBx* -af{y-x*) + ly'^By--h{y) | y e B(a;*, f) n c| , (3.12) 

for < f < r. Since h is continuous differentiable, the function inside the brackets is also continuous 
differentiable. Applying the first order necessary optimality condition for the smooth optimization 
problem (|3.12p (see, e.g., [H]), we obtain: 

Oe{B + B)x* -a- ^h{x*) + ^cnB(..,.-)(^*)- 



vcnB(..,.-)(^*) = Nc{x*) n ^S(,.,,)(x*) = Nc{x*). 



However, N^- 



Let d5c'{x) denote the subdifferential of the indicator function 6c of C at x. One has dSc{x) — 
Nc{x). Since matrix B is symmetric and positive definite, if we define gi{x) := ^x"^ Bx + Sc{x) then 
dgi{x) = Bx + d5c(x) and this mapping is maximal monotone. Consequently, := (/ + cdgi)~^ 
is well-defined and single valued, where / is the identity mapping (see I23j '). 

The following proposition provides a necessary and sufficient condition for a stationary point of 
(jncMVIPp . 

Proposition 3.2. A necessary and sufficient condition for a point x ^ C to be a stationary point 
to problem (jncMVIPp is: 

x = {I + cdgiy^ (x - c{Bx - a) + c\7h{x)j , (3.13) 
where c > and I stands for the identity mapping. 

Proof. Since gi is proper closed convex, the inverse (/ + dgi) ^ is single valued and defined every- 
where [22]. Thus X satisfies (j3.13p if and only if x — c{Bx ~ a) + cVh{x) ^ {I + cdgi){x). Moreover, 
since Nc{x) is a cone and dgi{x) — Bx + dSc{x) = Bx + Nc{x), the latter inclusion is equivalent 
to e Bx — a + Bx — Vh{x) + Nc{x), which shows that a; is a stationary point of (IncMVIPp . □ 

Now, if we define ydx) :— x ~ c{Bx ~ a) + c\7h{x) and 

Scix) {I + cdgiy\x - c{Bx - a) + cVhix)) , (3.14) 

then, it follows from Proposition l3 . 2l that x — Sc{x). Therefore, every stationary point x of (jncMVIPp 
is a fixed-point of S'c(-). To compute Sc{x), it requires to solve the following strongly convex quadratic 
problem over a convex set: 

mml^^y'^By+^J\y-y,{x)\\^ I x G c| , (3.15) 

This problem has a unique solution for any c > 0. 

Finally, we introduce the following concept, which will be used in the sequel. For a given tolerance 
e > 0, a point x* € C is said to be an e-stationary point to (jncMVIPp if 

D(P{x*;d) > -e, \/d e J"c(a;*), \\d\\ = 1. (3.16) 
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4 Gradient mapping and its properties 

By substituting ydx) into p.lSp . after a simple rearrangement, we can write problem p.lSp as 



min <{ ^y^By +[Bx~a~ \7h{x)f{y - a:) + J-||y - x\\^ \yeC 



(4.1) 



Now, we consider the following mappings: 

mc{x; y) ^-y^^ By+[Bx-&-yh{x)]^{y - x)-h{x) + ^\\y - xf, (4.2) 
Z zc 

and Sc{x) := argmin {mc(x; y) | y G C} . (4-3) 
Then, since problem (|4.3p is strongly convex, Sc{x) is well-defined and single- valued. Let us define 



Gcix) -[a; - Sc{x)]. 



(4.4) 



The mapping Gc(-) is referred as a gradient-type mapping of (j3.3l) |17] . Appyling the optimality 
condition for (14.31) we have 



Bsc{x)+Bx-a-Wh{x)-Gc{x) {y- Sc{x))>0, Vy E C. 



(4.5) 



From now on, we further suppose that the cost function h is Lipschitz continuous differentiable 
on C with a Lipschitz constant > 0, i.e. 



\\Wh{x)-Wh{y)\\<Lh\\x~y\\, Vx, y S C. 
By using the mean-valued theorem, it is easy to show that the condition (j4.6p implies 



\h{y)-h{x)-Vh{xfiy-x)\ < hh\\y~xf, yx,y e C. 

The following lemma shows some properties of D(f>{-; •). 
Lemma 4.1. For any x E C , we have 

D^s,{x);x-s,{x)) > ^~''^\'^™ \\GrXx)r, 
Z?0(s,(x);y-s,(a;))>-[l+c(i„ + ||S||)]||Ge(a;)||||2/-s,(x)||, Vj/eC. 

As a consequence, for any d e J-cisc{x)) with \\d\\ — 1, we have 

D4>{s,{x)-d) > -[l + ciLh + \\B\\)]\\Gcix)\\. 

Proof. From the definition of D(j) in p.lOp . we have 

r ~ 1 ^ 

D(f)(sc{x); X — Sc{x))= Bsc{x)~a+Bsc{x)~\'h{sc{x)) {x — Sc{x)) 

1 T 

Bx — (5 — \7h{x) + Bsc{x) (x — Sc{x)) 

-1 T 

Bsc{x) — a — \7h{sc{x)) — Bx — a + \7h{x) {sc{x) — x) 



(4.6) 



(4.7) 



(4.8) 
(4.9) 

(4.10) 



(4.11) 



> 



Bx-a- Wh{x) + Bsc{x) {x ~ Sc{x)) - {Lh + \\B\\)\\x ~ Sc(a;)|p 



Substituting (|4.5p into (|4.11l) we obtain 



Dcf>{sc{x);x-Sc{x)) > {- - [Lh + \\B\\)]\\x ~ s,{x)\\' 



1-c{Lh + \\B\\). 



Gcix)\\ 
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which proves 

Using again (|4.5p and (|4.6p we have 



D(j){sc{x);y - Sc{x)) ^ Bsc{x) - a + Bsc{x) - Vh{scix) {y - sdx)) 



Bsc{x) - a- Vh{sc{x)) (y ~ Sc{x)) + {Bsc{x)Y{y - x) 



> 



Bsc{x) ~ a -Vh{sc{x)) {y - Sc{x)) 



Bx — \'h{x) H — {sc{x) — x) 



(scix) - y) 



= Bsc{x)-Vh{sc{x))-{Bx-a)+Wh{x) {y - Sc{x)) + Gc{xf {sc{x) - y) 
> ^{Lh + \\B\\)\\x - s,{x)\\ \\y - s,{x)\\ ~ \\G,{x)\\ \\y - s,{x)\\ 



> - 



l + c(L„ + ||B||) \\G,{x)\\\\y ^ s,{x)\\, 



which proves (14.91) . 

By the convexity of C, there exists t >Q such that Sc{x) + td E C, where \\d\\ = 1. If we substitute 
y := Sc{x) + td E C into (14. 9p then we get 



D<i)(x-td) > -t{l + cLfh + c\\B\\)\\Gc{x)\\. 



(4.12) 



If i = then (j4.9p automatically holds. If t > then by the linearity of D with respect to the second 
argument, we divide both sides of (|4.12p by i > to get (|49| . □ 

Remark 4.1. For a fixed x G C, if we define edx) :— \\Gc{x)\\ and rdx) := \\x — Sc(a;)|l then edx) 
decreases in c and rdx) increases in c. 

Indeed, let q{y, c) := {Bx - of {y - x) + ^y'^By - h{x) ~ \/h{x)'^{y ~ x) + ^\\y - x\\^. Then 
q is convex jointly in two arguments y and c. Thus uj{c) :— miuj^gc g(y, c) is convex. It is easy to 
see that w'(c) = — i||G'c(2;)|P increases in c. Hence, edx) decreases in c. If we replace c by 1/c in 
q{y,c), this function becomes concave in c, then by the same argument as cij(-), we conclude that 
rdx) increases in c. 

Since B and B are symmetric, we consider a potential function defined as follows: 



7(2;) :— —x^Bx + -x'^ Bx — X — h{x), 



(4.13) 



Then, 7 is nonconvex but Lipschitz continuous differentiable. We have the following statement. 
Lemma 4.2. For x,y G G, we have 



mdx;sdx))+x Bx - a x < j{x) - -\\Gdx)\\' 



Moreover, if c{Lh + ||-B||) < 1 then 



mdx; sdx)) + -X Bx - a x >j{sdx)). 



(4.14) 



(4.15) 



Proof. It is obvious from the definition of mc(x;x) that (l){x) — mc{x;x) + ^x Bx — a x. Since 
mdx; ■) is strongly convex quadratic with modulus ^, using (j4.13p we have 



7(2;) — mdx; sdx)) — -x^ Bx + dd'x = mdx; x) — mdx; sdx)) 

>h^\x-sdx)r^^-\\Gdx)\\\ 



8 



which proves (|4.14p . 

To prove (|4.15l) . from (|4.7p and the definition of 7 we have 

j{sc{x)) — mc{x; Sc{x)) — -x'^ Bx + oF x — — Sc{x)Bsc{x) — x^Bx ~ 2{Bx)'^ {sc{x) — x) 

- h{sc{x)) + h{x) + Vhixfiscix) - x) - ^\\sc{x) - xf 

< - xfBMx) <fil|-:il||s,(i) - if 

<- '-'"^;""^"' llM.)-.f- ° 

By assumption c{Lh + \\B\\) < 1, we obtain (j4.15p . □ 
If we combine the inequalities (j4.15p and (14.141) in Lemma [4.21 then: 

7(5c(x))<7(x)-^||G,(x)|p. (4.16) 

This inequahty plays an important role in proving the convergence of the splitting proximal point 
algorithm in the section. 

For a given starting point e C, let us define the level set of 7 with respect to C as 

C^ijix")):^{xeC\j{x)<jix")}. (4.17) 

From it is obvious that if a;° e £^(7(2;°)) then Sc{x") G £^(7(2;°)) provided that c(Lh + ||_B||) < 

1. 

5 A splitting proximal point algorithm and its convergence 



Proposition 13.21 suggests that a proximal point method can be applied to find a stationary point of 
(lncMVIP[) . For the implementation purpose, the proximal mapping defined by p.l4p is extracted 
to the expression p. 151) . The splitting proximal point algorithm constructs an iterative sequence as 
follows: 

Algorithm 1. (The splitting proximal algorithm) 

Initialization: Choose a positive number cq > 0. Find an initial point x''' E C and set k := 0. 
Iteration k: For a given x'^, execute the three steps below. 

Step 1: Evaluate Vh{x'^) and set yk := x*-' — Ck{Bx'^ — a) + CkVh{x^). 

Step 2: Compute x^^^ by solving the following convex quadratic program over a convex set: 

min|ix^Ba; + ^|ly-/f I yecj. (5.1) 

Step 3: If ||a;'^+-'^ — x^\\ < e for a given tolerance e > then terminate, x'^ is an e-stationary 
point of (jncMVIP|) . Otherwise, update ct and increase fe by 1 and go back to Step 1. 

In Algorithm [T] we left unspecified the way to update c^. If the Lipschitz constant Lh is provided 
then we can choose — for all fc, where :— Lh + ||^||- Otherwise, a line-search procedure 
can be used to update Ck- The latter procedure is briefly described as follows. First, we choose two 
constants c and c such that c > and < c < +00. Then we perform the following steps. 

• Given a constant Tc G (0, 1). Choose an initial value of c in [c, c]. 
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• Compute Scix'^)- While the decreasing condition 



liscix"}) < me(a;^ + ^{x'^fBx'' - a'^x\ (5.2) 



does not satisly, increase c by c := TcC and recompute Sc{x''). 
• Set Cfc+i := c. 



Now, we define Aa;*^ :— x''^-^ — x'^ and 



4 mm — . (5.3) 

o<i<fc 2cj 



The convergence of the splitting proximal point algorithm is stated as foUows. 

Theorem 5.1. Suppose that the function h is Lipschitz continuous dijferentiable on C with a Lip- 
schitz constant Lh > 0. Suppose further that for a given a;" G C the level set £^(7(0;°)) is bounded 
(particularly, C is bounded). Then the sequence {x'^}k>Q generated by Algorithmic starting from x^ 
satisfies: 

(7(^0) -7) 

5k < ' \ ' Vfc > 0, (5.4) 

A; + 1 

where 7 := inf 7(2^)- Moreover, for any d G J-c(x^'') with \\d\\ ~ 1, we have 



- 2(7(a;0)-7) 

Dcj^ix'^ ; d) > -(1 + cU + cIlSlI) W (5.5) 

where ik is the index such that Qj. ||Gcj^ (a;*'' )||^ = A^. 

As a consequence, if the sequence {x''} generated by (j5.ip is bounded, then every limit point of 
this sequence is a stationary point of (jncMVIPp . The set of limit points is connected and if it is 
finite then the whole sequence {x^} converges to a stationary point of (|ncMVlP|) . 

Proof. Since £^(7(0:'^)) is bounded by assumption, we have 7 :~ uif^eCyi-yix")) is well-defined 
due to the continuity of 7 and the closedness and nonemptiness of >C^(7(x°)) (since x^ G £^(7(x*')). 
From Step 3 of Algorithm [Tl if either the constant parameter Ck — -r— or the line search procedure 
is used then it implies 

^(2;fe+i) + — llx'^'+i - x'^f < m,,(x''-+i) + {Bx'' - afx'' < -fix''), Vfc > 0. (5.6) 
2ck 

Note that the whole sequence {x'^} is contained in Cj{'j{x'^)). Rearrange and sum up these inequal- 
ities for fc = to fc = ii' we get 

^' 1 

k—Q 

Then the inequality (j5.4l) directly follows from the definition of S in (j5.3l) . Combining (15.41) and 
(|4.10p in Lemma liH we obtain (|5.5p . 

To prove the remainder, taking into account Remark 14.11 and then passing to the limit as k tends 
to 00 the resulting inequality of (|5.7p . we get 



y -\\x^+' -x'^w^ <+oo. 



k=0 

Since c < -l-oo, this inequality implies that limfc^oo ||a;'^+^ — x'^W = 0. Therefore, the set of limit 
points is connected. Combine this relation and the assumptions of boundedness of {x''} it is easy 
to show that every limit point of {x'^} is a stationary point of (jncMVIPp . When the set of the limit 
points is finite, the last statement of the theorem is proved similarly using the same technique as in 
[20,[Chapt. 28]. □ 
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Remark 5.1. For a given tolerance e > 0, according to Theorem 15.11 the number of iterations k 
to get an e-stationary point is O(e^). Consequently, the worst-case complexity of Algorithm [T] is 
Oil/Vk). 

6 Numerical test 



In this section, we consider to numerical examples involving concave cost functions. The aim of 
these examples is to estimate the number of iterations of Algorithm [1] in a certain case compared 
to the worst-case complexity given in Theorem 15. II In addition, we also test the time profile of the 
algorithm when the size of problem increases. 

The algorithm is implemented in Matlab 7.8.0 (R2009a) running on a Pentium IV PC desktop 
with 2.6GHz and 512Mb RAM. We assume that the feasible set C of (|ncMVIP|) is a box in E". 
Therefore, the convex problem p.l5p reduces to quadratic programming. We solve this problem by 
using the quadprog solver (a built-in Matlab solver). 

Example 1. Suppose that the cost function hi{xi) of the firm i is given as hi(xi) = c'^ + Oi ln(l-f r^Xi), 
where > is the ceiling cost, Ci > and > are given. The function h becomes 



h{x) = c" + ^ c,; ln(l + nxi) = c° + In J|(l + nx^Y 



(6.1) 



where c° = J27=i '^i - '^^ obvious that hi is well-defined if a;^ > and h[{xi) 
that h is differentiable on C = W] and 



which implies 



Vh{x) = ( 



ciri 



1 + rixi ' 



1 -t- TyiXjl 



(6.2) 



Since h'l{xi) — —Cirf/{l+riXi)'^, we have |ft,"(a;i)| < cirf and h is concave. Moreover, Vft. is Lipschitz 
continuous with the Lipschitz constant :— maxjcirf | i = 1, . . . , n}. 

In this example, we choose /3 = 0.1 > 0, a = 10, = 2, Ci = 1.5 for all i = l,...,n, and 
rt = 1 + uji, where uJi is randomly generated in (0,1) {i — 1, . . . , n). The strategy set of the firm i is 
defined by Ci := [0, 10] for alH = 1, . . . , n. 

We test Algorithm [T] for problem (jncMVIP[) with the size increasing from 10 to 1000. The 
tolerance e is 10"'^. The number of iterations as well as the CPU time with respect to the size of 
problem is visualized in Figl. and Fig2., respectively. 





Figl. Number of iterations depending on n [Ex. 1] Fig2. CPU time depending on n [Ex. 1[ 

From (15.51) of Theorem 15. 1[ it implies that the number of iterations k to reach an e-stationary 
point depends on the structure of the function h and the value 7(a;°) — 7, Lh and \\B\\. Since ft, is a 
logarithm function, the value of h slowly increases in n, while the Lipschitz constant L/i < 1.5 x 2^ — 6 
for all n and the norm = (n — l)/3. Consequently, the worst-case complexity bound increases 
almost linearly in n. As can be seen from the first figure, the number of iterations increases with 
a small slope when the size of problem grows up. The curvature of this graph stays below a linear 
line generated by the worst-case complexity bound. The CPU time also increases almost linearly in 
the size of problem. 
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Example 2. In this example, we choose the cost function hi as hi{xi) = c° — c^e"'"'^', where 
c° > Ci > and > given. It is easy to see that h'-{xi) = —Ciafe'^'^'- < 0, then hi is concave. 
Since hi is differentiable on M, it means that h is differentiable on and V/i is expressed by 

Vh{x) = (cirie-'^i^\.-- ,c„r„e-'-"^")^. (6.3) 

We have \hi{xi)\ < Cirf for alH = 1, . . . , n, thus V/i is Lipschitz continuous on R" with the Lipschitz 
constant Lh := maxjcir^ | 1 < * < n}. 

To compare with the previous example, we choose the value of the parameters a and /3 as in 
Example 1. The parameters and Ci are given by = 4 and Ci = 2 for all i = 1, . . . ,n. The 
parameter 0.1 + O.lrandi, where rand^ is generated randomly in (0, 1). 

1400 p 
1200 - 
1000 - 
I 800 - 
600 - 
400 - 
200 - 
o ^ 

100 200 300 400 500 600 700 800 900 1000 100 200 300 400 500 600 700 800 900 1000 

n [size ol problem] n [size o! problem] 

Fig3. Iterations depending on n [Ex. 2] Fig4. CPU time depending on n [Ex. 2] 

We also test Algorithm [T] for the problem size from 10 to 1000. The number of iterations and 
the CPU time are plotted in Fig3. and Fig4., respectively. Since the function 7 rapidly increases 
in n compared to the previous case, the number of iteration also increases. Consequently, the CPU 
time respectively increases. 
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